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Abstract 

It  Is  well  known  that  explosive-based  propellants  are  susceptible  to 
detonation  from  the  controlled  deflagration  mode  of  combustion.—’ ^fn^certain 
instances  (i.e.  when  the  propellant/explosive  is  fragmented)  the  likelihood  of 
a  catastrophic  event  is  greater.  Fragment  size,  gas  permeability  through  the 
packed  bed,  chemical  decomposition  rate  and  product  gas  confinement  all  play 
an  important  role  in  determining  whether  a  convective  deflagration  wave  will 

CT" 

make  a  transition  to  a  steady  state  detonation  wave.'^In  some  instances  a 
confined  zone  of  granulated  propellant  adjacent  to  a  zone  of  cast  propellant 
can  provide  a  rapid  enough  pressure-rise  rate  to  shock  initiate  the  cast 
material.  If  the  cast  propellant  has  voids,  the  detonation  will  initiate  at 
some  location  ahead  of  the  granulated  bed/cast  material  interface.-v 


r — 

^-^This  report  is  a  summary  of  the  research  activities  that  focus  on  the 
analysis  and  modeling  of  the  physics  of  such  highly  transient  flows^The 
report  includes  (as  Appendices)  two  technical  papers  prepared  during  the  past 
year  on  the  topic  of  DSOT  (Deflagration  to  Shock  to  Detonation  Transition). 
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1.  Introduction 

The  probability  of  a  detonation  occurring  in  a  solid  propellant  rocket 
motor  greatly  Increases  when  secondary  high-explosives  (HE)  are  used  as 
constituents  in  the  propellant  mixture.  Octogen  (HMX)  and  cyclotrimethylene 
trlnitramlne  (RDX)  are  the  most  commonly  used  nitramine  HE  in  todays 
propellant  formulations.  By  Increasing  the  explosive  content  of  solid 
propellants  the  overall  performance  of  the  system  increases  as  well  as  the 
sensitivity  to  detonation  by  shock  initiation.  Since  the  early  1970's 
laboratories  at  both  the  Air  Force  and  Navy  have  researched  and  tested  solid 
propellant  formulations  containing  secondary  explosives  such  a  RDX  and  HHX. 
As  will  be  discussed  briefly,  a  detonation  hazard  exists  when  burning  these 
high-energy  propellants,  which  is  not  present  during  the  combustion  of 
"standard"  composite  and  composite-double  base  propellants. 

By  burning  more  energetic  propellant  constituents  the  specific  impulse 
and  thus  overall  performance  of  a  solid  rocket  motor  can  be  increased. 
However,  one  disadvantage  which  comes  with  using  these  explosive-based 
propellant  mixtures  is  the  hazard  of  Deflagration  to  Shock  to  Detonation 
Transition  (DSDT).  A  DSDT  event  is  when  a  controlled  subsonic  deflagration 
wave  makes  a  transition  to  a  high  order  detonation.  The  result  is  total 
destruction  of  the  motor  assembly. 

In  the  literature  (reviewed  in  Appendices  A  and  B)  this  process  is  more 
commonly  referred  to  as  Deflagration  to  Detonation  Transition  (DDT).  Our 
research  group  at  Illinois  coined  the  description  DSDT  in  order  to  emphasize 
that  Shock  to  Detonation  Transition  (SDT)  is  the  final  step  in  the  process. 
Both  terms,  DDT  and  DSDT,  will  be  used  interchangeably  throughout  the  text. 


DSOT  in  Granulated  Solid  Propellant 
In  the  papers  Included  as  Appendices  A  and  B,  three  different  flow 
processes  have  been  Identified  as  DDT  (DSDT).  Two  are  real  possibilities  In 
the  rocket  motor  environment  and  the  third  Is  only  possible  when  using  a  high 
order  cast  explosive  under  extreme  confinement.  (It  is  also  unlikely  that  a 
solid  rocket  motor  would  be  cast  with  pure  explosive.)  But  studying  this 
third  scenerio  provides  basic  information  on  detonation  initiation 
mechanisms.  This  section  has  been  included  in  order  to  differentiate  between 
the  three  cases  and  provide  the  reader  with  a  basic  understanding  of  the 
sequence  of  event  leading  to  detonation  in  each  of  them. 

Case  1 

What  will  be  referred  to  as  "DDT-Case  1"  involves  a  transition  to 
detonation  occurring  within  a  bed  of  totally  porous,  granulated  propellant. 
Consider  the  rocket  motor  shown  in  Fig.  I.  For  illustration  purposes,  a 
center-burning  configuration  is  shown.  It  is  hypothesized  that  the  normal 
burning  process  of  the  cast  solid  propellant  in  the  rocket  motor  is  disturbed 
by  an  abnormality  such  as  a  crack  (see  Fig.  1)  in  the  propellant  grain  [1,2], 
thus  providing  the  granulated  region.  Figure  2  shows  an  enlargement  of  the 
granulated  region.  The  fracture  could  be  the  result  of  a  handling  accident 
during  shipment  or  nozzle  failure  during  operation.  Because  of  the  increased 
surface-to-volume  ratio  of  the  resulting  propellant  fragments,  product  gas 
generation  increases  beyond  the  level  necessary  for  steady-state  motor 
operation.  That  is,  the  mass  produced  in  the  chamber  exceeds  the  mass  leaving 
the  nozzle.  Pressure  gradients,  developed  as  a  result  of  localized  burning, 
drive  the  hot  product  gases  into  the  cracks  developed  in  the  upstream 
propellant.  As  a  consequence  of  this  unsteady  flow  process,  convective  heat 


DSDT  HAZARD  IN  ROCKET  MOTOR 


Fig.  1 


Sketch  of  solid  propellant 
rocket  motor  with  cracked 
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transfer  from  the  hot  gas  to  the  unreacted  propellant  will  ignite  more  of  the 
upstream  propellant  particles.  Under  certain  circumstances  the  process  can  be 
self-propagating.  That  is,  as  the  propellant  decomposes,  the  pressure  wave 
strengthens,  which  in  time  leads  to  the  ignition  of  more  propellant.  This 
accelerated  convective  burning  can  eventually  lead  to  shock  compression  of  the 
upstream  propellant  and  a  possible  detonation  transition. 

Figure  3  illustrates  the  general  form  of  a  pressure-distance  profile  once 
the  steady  state  detonation  solution  is  obtained.  For  illustration  purposes, 
the  profile  is  superimposed  on  a  sketch  of  a  granulated  bed.  The  shock  front 
is  followed  by  a  narrow  reaction  zone  (not  to  scale  in  Fig.  3)  which  is 
followed  by  an  expansion  zone  consisting  of  all  product  gases.  Appendix  A, 
“Analysis  of  Deflagration  to  Shock  to  Detonation  Transition  (DSDT)  in  Porous 
Energetic  Solid  Propellants,"  by  Butler  and  Krier,  is  a  detailed  study  of  the 
"DDT-Case  1"  scenerio.  The  reader  should  consult  Appendix  A  for  a  discussion 
on  the  governing  equations  and  predicted  DSDT  results. 

Case  2 

A  second  DDT  scenerio,  "DDT-Case  2,"  involves  the  rapid  local 
pressurization  in  a  region  of  granulated  propellant  as  the  necessary  impetus 
to  shock  initiate  an  adjacent  region  of  cast  but  void-containing  propellant. 
This  is  illustrated  in  Fig.  4.  The  cast  material  (Zone  1)  can  contain  'blind' 
pores,  but  is  assumed  to  be  impermeable  to  the  flow  of  hot  gases  form  the 
granular  zone  (Zone  2).  This  implies  that,  unlike  the  first  DDT  scenerio 
discussed,  only  stress  waves  can  be  transmitted  upstream  of  the  reaction 
zone.  A  second  characteristic  of  Case  2  is  that  the  length  of  the  granular 
bed  is  less  than  the  detonation  run-up  length,  The  important  point  here 

is  that  although  the  granular  bed  is  shorter  than  the  critical  detonation  run- 
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-Expansion  Zone - H^^Zone*'00 


g.  3  Schematic  of  Packed  Bed  after  Transition  to  Detonation. 
Ignition  region  has  collapsed  to  a  thin  zone  and  is 
followed  by  an  all  gas  expansion  zone. 


Zone  2  Zone  1 


Burning  k 

Solid 


Unreacted  Solid 


x=L 


Fig.  4  Sketch  of  granulated  bed/cast  bed 
configuration  (Case  2)  for  DSDT. 
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up  length,  the  impermeable  material  can  still  detonate.  DDT-Case  2  was 
studied  In  detail  by  Cudak,  Krier  and  Butler  [3]  which  is  attached  as 
Appendix  B. 

Figure  5  shows  the  pressure-time  profiles  acting  on  the  granular  bed/cast 
solid  Interface  for  various  propellant  sizes  [3].  As  shown  in  the  figure,  the 
pressure-rise  rate  is  strongly  dependent  upon  the  grain  size.  The  larger  the 
particle  size,  the  slower  the  pressure-rise  rate.  The  cases  studied  in  Ref.  3 
had  rise-time  of  0<t*<30ys.  This  pressure-rise  function  is  used  as  a  boundary 
condition  for  the  unsteady  flow  analysis  upstream  of  the  interface. 

The  governing  equations  defining  stress  wave  propagation  and  subsequent 
detonation  upstream  of  the  interface  are  the  usual  flow  equations  for  one¬ 
dimensional  flow  in  a  continuous  medium.  In  a  Lagrangian  formulation,  one 


writes. 


vt  =  vux 


ut  =  -vPx 


e*  *  -Pvt  +  Q 


where  v  is  specific  volume,  u  particle  velocity,  P  pressure,  e  specific 
internal  energy,  and  Q  chemical  energy  release  rate.  The  subscripts  ' x '  and 
' t '  represent  partial  derivatives  with  respect  to  space  and  time 
respectively.  Here,  one  is  conserving  the  mixture  properties,  i.e.  both  solid 
and  gas. 

In  addition  to  the  conservation  equations  (Eqs.  1-3),  an  equation  of 
state,  P  -  P(e,v),  (Eq.4),  and  the  reciprocal  caloric  equation  of  state 
e  *  e(v,T),  (Eq.  5),  are  needed  for  both  the  solid  and  product  gas  phases. 
Finally,  a  stress-void  volume  pore  collapse  relation 


Pressure  (dynes /cm2) 


and  energy  partition  equation 


ET  =  Es  +  Eg  (7) 

are  required  to  complete  the  mathematical  description  of  the  flow  behavior 
prior  to  reaction.  In  Eq.  7,  E  represents  total  energy,  subscript  's'  the 
solid  phase  and  subscript  ' g*  the  gas  phase.  Appendix  B  provides  additional 
information  concerning  these  equations. 

Once  reaction  commences,  a  chemical  decomposition  rate 

-  f(W,E*,T*)  (8) 

is  also  needed  to  define  the  complete  problem.  Here,  W  represents  an 
undecomposed  mass  fraction,  E*  an  activation  energy  and  T*  represents  a 
reaction  temperature. 

Experimental  work  [4]  has  shown  that  material  samples  containing  voids 
and  density  irregularities  will  undergo  shock  to  detonation  transitions  at 
much  lower  shock  pressures  than  a  homogeneous  sample  of  the  same  material. 
Shock  waves  not  strong  enough  to  raise  the  bulk  temperature  of  the  material 
above  the  thermal  explosion  level  are  distorted  by  the  density  discontinuities 
of  the  porous  material  and  subsequently  superheat  the  material  in  these 
localized  regions  above  the  explosion  level.  Following  this,  the  material 
reacts  and  strenghtens  the  leading  shock  wave  which  causes  the  transition  to  a 
detonation.  These  "hot  spots",  as  they  are  referred  to,  are  an  initiating 
mechanism  in  porous  reactive  material  [5], 
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The  emphasis  of  the  work  cited  in  Appendix  B  was  to  show  that  ramp¬ 
loading  porous  materials  to  relatively  low  peak  stresses  can  lead  to 
detonation.  In  most  cases,  the  bed  detonated  at  a  location  downstream  of  the 
stress  loaded  boundary  and  an  energy  sustaining  retonation  wave  was  predicted 
to  travel  back  towards  the  Interface. 

Figure  6  presents  the  pressure  (stress  wave)  build-up  history  in  a 
reactive  solid  with  blind  voids  (a  =  1.117),  being  pressurized  to  2  GPA  at  a 
linear  rate,  taking  10  psec. 

In  each  portion  of  Figure  6  [parts  (a)  -  (j)]  we  show  the  predicted 
pressure  (stress)  at  the  fixed  time  (shown  in  the  upper  left  hand  corner  of 
the  figures).  Note,  even  as  early  as  8  psec  the  material  has  been  compressed 
by  several  millimeters,  as  evidenced  by  the  'gap'  in  the  p  -  x  variation  near 
x  =  0.  At  t  =  16  psec  a  steepend  shock  is  clearly  seen  at  x  =  3.4  cm.  At 
18  psec  a  detonation  wave  is  predicted  to  begin  traveling  both  forward  and 
rearward. 

Within  another  microsecond  these  oppositely-traveling  reactive  shocks  are 
almost  several  centimeters  apart,  since  the  predicted  detonation  and 
retonation  waves  move  at  speeds  of  8.3  and  8.7  mm/psec,  respectively.  The 
retonation  reflects  from  the  front  wall  (interface)  and  a  more  complicated 
pressure  distribution  is  shown  in  Figs.  6(h),  6(i)  and  6(j). 

The  calculations  shown  in  Figs.  6(a)  -  6(j)  considered  an  explosive  (HMX) 
with  an  energy  of  6.84  Mjoule/kg  having  a  porosity  of  10.5%  (a  =  1.1176).  The 
material  of  this  density,  according  to  thermodynamic-equilibrium  calculations 
is  predicted  to  detonate  at  a  C-J  pressure  of  27.9  GPa.  Note  that  our 
transient  flow  calculations  predicted  an  eventual  steady  C-J  pressure  for  the 
detonation  to  be  27  GPa,  a  very  good  agreement. 


Pressure  (GPa)  Pressure  (GPa) 


30. 


20. 


10. 


f-  8.0 
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Fig.  6  Pressure  (stress)  distribution  at 

discrete  times;  time  event  encompasses 
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Figure  7  presents  the  locus  of  the  stress  wave  (x-t  diagram),  shown  as 
the  dotted  curve.  At  17.8  ysec,  at  a  location  of  2.43  cm,  a  detonation  (and 
retonatlon)  Is  predicted.  The  slope  In  the  x  -  t  plane  Is,  of  course,  the 
reaction  front  speed,  and  It  Is  predicted  to  be  8.3  mm/ysec  for  the 
detonation,  and  8.8  mm/ysec  for  the  returning  detonation  (retonatlon),  which 
travels  into  a  p re- compressed  (higher  density)  explosive. 

Appendix  B  shows  more  results  and  provides  details  of  the  analysis.  It 
Is  clear  that  our  analysis  is  pertinent  to  the  detonation  physics  responsible 
for  DSDT  is  porous  explosives.  Of  course,  additional  work  Is  In  order,  since 
the  sensitivity  of  the  many  assumptions  on  the  predicted  detonation  run-up 
length  has  yet  to  be  determined.  For  example,  the  variation  In  chemical 
activation  energy  during  the  rapid  reaction  process  that  supports  the 
detonation  can  have  a  significant  effect  in  the  predicted  run-up  lengths. 

Case  3 

In  both  DSDT  processes  described  above,  the  transition  to  detonation  was 
the  result  of  convective  burning  in  a  packed  bed  of  reactive  particles.  In 
Case  1  the  detonation  transition  occurred  within  the  granulated  bed  and  in 
Case  2  the  detonation  occurred  in  a  cast,  porous  material  upstream  of  the 
granulated  bed.  Nevertheless,  in  both  cases  the  rapid  pressurization  rate  due 
to  the  reacting  fragmented  bed  provided  the  driving  force  necessary  to  shock 
initiate  the  material  upstream.  The  third  type  of  DDT  discussed  in  the 
literature,  "DDT-Case  3,"  results  from  end-burning  (conductive  combustion)  a 
confined  cast  explosive  (Fig.  8)  which  is  impermeable  to  the  flow  of  hot 
product  gases.  Although  the  manner  in  which  the  deflagration  wave  traverses 
the  explosives  is  different  from  the  first  two  cases,  the  end  result  (steady 
detonation  wave)  is  the  same. 
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Fig.  7  Locus  of  stress  front  (dashed  line)  and 

the  detonation  and  retonation  front  (solid  line) 
for  pressure  history  shown  in  Figures  6(a)-6(j). 


Figure  8  is  an  illustration  of  the  test  configuration  used  in  this  type 
of  DOT  experiment.  Here,  the  region  labeled  'Zone  1'  is  a  cast,  voidless 
explosive  and  'Zone  2'  is  occupied  by  the  product  gases  generated  in  the 
propellant  combustion.  The  sequence  of  events  for  this  DDT  begins  with  the 
thermal  ignition  of  the  explosive  at  the  location  x=0.  This  is  then  followed 
by  the  pressurization  of  the  gas  volume,  (Zone  2),  increased  regression  rate 
(dx/dt)  of  the  burning  material,  stressing  of  the  unreacted  solid,  shock 
formation  ahead  of  the  ignition  front  and  eventual  transition  to  detonation  at 
the  location  of  the  upstream  shock  front  formation.  The  driving  force 
provided  by  the  confined  gases  in  Zone  2  is  analogous  to  a  moving  piston  with 
a  prescribed  velocity-time  profile  acting  on  the  x=0  boundary,  continually 
increasing  the  stress  level  in  Zone  1.  Although  Case  3  has  been  shown  to  be  a 
viable  method  for  obtaining  DDT  in  cast  explosives,  it  will  not  cause  DDT  in 
cast  propellants.  Very  simply,  conductive  combustion  alone  will  not  provide  a 
rapid  enough  pressure-rise  rate  to  shock  initiate  cast  propellant. 

1.3  Conclusions 

In  the  work  presented  here  and  in  our  papers  attached  as  Appendices,  we 
have  identified  three  different  ways  in  which  DSDT  can  occur  in  a  solid  rocket 
motor  which  has  been  damaged.  The  first  two  (Cases  1  and  2)  are  real 
possibilities  in  actual  motors  and  the  third  is  highly  unlikely.  Case  3 
requires  a  perfect,  voidless  high  energy  reactive  solid,  in  order  to  undergo  a 
transition  to  detonation  to  the  bulk  material  heating  by  rapid  shock 
compression. 

In  Appendix  A  we  showed  that  DDT  can  occur  when  a  confined,  but  granular 
bed  of  high-energy,  high  reaction  rate  material  is  ignited  at  one  end.  The 
transient  build  up  to  detonation  was  only  possible  because  of  the  convective 
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heat  transfer  mechanism  from  the  hot  gas  to  the  unreacted  propellant,  early  in 

the  process.  This  DDT  process  was  identified  as  DDT-Case  1.  ; _ 

DDT-Case  2,  highlighted  in  Appendix  B,  occurs  when  the  pressure-rise  rate 
from  a  bed  of  granular  propellant  drives  an  adjacent  porous  cast  material  to 

detonation.  As  shown  in  Appendix  B,  this  can  occur  at  relatively  low  stress  . _ . 

levels  in  the  cast  material. 
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SUMMARY 

It  is  well  known  that  expl os  1 ve-hased  propellants  are  susceptible  to  detonation 
from  the  controlled  deflagration  mode  of  combustion.  In  certain  instances  (i.e.  when 
the  propel  1  an t / e* pi  os i ve  is  fragmented)  the  likelihood  of  a  catastrophic  event  is 
greater.  Fragment  size,  gas  permeability  through  the  packed  bed,  chemical  decompo¬ 
sition  rate  and  product  gas  confinement  all  play  an  important  role  in  determining 
whether  a  convective  deflagration  wave  will  make  a  transition  to  a  steady  state  deto¬ 
nation  wave.  The  work  presented  here  represents  an  effort  to  analyze  and  numerically 
model  the  governing  equations  defining  such  a  transient  two-phase  reactive  flow  pro¬ 
cess.  Computer  generated  results  will  be  presented  showing  the  development  of  a 
steady  state  detonation  wave  from  an  accelerating  deflagration  wave. 
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INTRODUCTION 

By  burning  more  energetic  propellant  constituents  such  as  the  secondary  explo¬ 
sives  cycl ic-trimethylenetrini ntremene  (RDX)  and  octogen  (hmx)  the  specific  impulse 
and  thus  overall  performance  of  a  solid  rocket  motor  can  be  greatly  increased.  How¬ 
ever,  one  disadvantage  which  comes  with  using  these  explosive-based  propellant  mix¬ 
tures  is  the  hazard  of  Oeflagration  to  Shock  to  Detonation  Transition  (OSOT) .  A  DSDT 
event  is  when  a  controlled  subsonic  deflagration  wave  maxes  a  transition  to  a  high- 
order  detonation.  The  result  is  usually  total  destruction  of  the  motor  assembly.  In 
a  rocket  motor  the  OSOT  is  usually  associated  with  propellant  fragmentation  prior  to 
the  event.  Although  an  actual  motor  has  never  been  instrumented  in  fine  enough 
detail  to  describe  every  event  leading  up  to  a  DSDT,  researchers  have  hypothesized 
what  the  sequence  of  prc-OSOT  events  might  he. 


'/////////////////////////////////^^ 

VI  aii  _ r-i _ • I  I  rv _ :  —  _  _  L  i  i  ■ !  a  _  r  J 


1 


AIL  Gas  Region 


Region  of  Unignited 
Particles 


! 


^77777777777777>77^77^ 


x=o 


Ignition  Region- 


X  =  L 


Fig.  1  Schematic  of  Packed  8ed  after  Transition  to  Detonation.  Ignition  Region  has 
Collapsed  to  a  Thin  Zone  and  is  followed  by  an  all  Gas  Region. 


As  stated  above,  the  first  step  in  the  OSOT  process  appears  to  be  fragmentation 
of  the  grain.  This  can  be  the  result  of  a  handling  accident  during  shipment  or  a 
nozzle  failure  during  operation.  Recause  of  the  increased  su rf ace- t o- vo 1 ume  ratio  of 
the  resulting  fragments,  product  gas  generation  increases  beyond  the  level  necessary 
for  steady-state  motor  operation.  That  is,  the  mass  produced  in  the  chamber  exceeds 
the  mass  leaving  the  nozzle.  Pressure  gradients,  developed  as  a  result  of  the  loca¬ 
lized  hurning,  drive  the  hot  product  gases  into  the  cracks  developed  in  the  upstream 
propellant.  As  a  consequence  of  this  flow  process,  convective  heat  transfer  from  the 
hot  gas  to  the  unreacted  propellant  will  ignite  more  of  the  grain  particles.  The 
accelerated  convective  burning  can  eventually  lead  to  shock  compression  of  the  up¬ 
stream  propellant  and  eventual  detonation.  Figure  1  is  an  illustration  of  a  packed 
bed  once  the  transition  has  occurred. 

Due  to  the  obvious  hazard  and  complexity  of  the  detonation  transition  process, 
there  are  limited  experimental  data  available  to  verify  the  proposed  model  describing 
this  dynamic  event.  It  is  nevertheless  appropriate  at  this  point  to  discuss  and 
reference  some  of  the  experimental  work  on  OSDT  prior  to  the  discussion  of  our 
modeling  effort. 

In  order  to  simulate  the  de f 1  a g ra t i on - 1 o- s hoc k - 1 o- de t on  a t i on  (OSDT)  event 
occurring  in  a  solid  rocket  motor  without  the  destruction  of  'u 1 1  scale  motor  assem¬ 
blies,  several  research  groups  [  l  -  3  1  have  studied  the  detonation  transition  present 
in  a  packed  bed  (10-20  cm  in  length)  of  granulated  p rope  1 1  a n t / e < o 1  os i ve .  The  deto¬ 
nation  tube  is  dosed  at  both  ends  and  instrumented  throughout  to  record  stress 
levels  and  ignition  front  loci  as  the  flame  propagates  from  the  igniter  end,  in  some 
cases  transitioning  to  a  steady  state  detonation  wave. 

By  analyzing  the  stress-time  and  ignition  front  data  obtained  in  these  types  of 
experiments,  the  authors  of  Reference  1  have  presented  a  scenerio  of  the  events 
leading  to  the  detonation  transition.  in  brief,  this  includes  the  propagation  of  a 
weak  compaction  wave  through  the  unreacted  explosive  due  to  the  igniter  blast,  con¬ 
vective  burning  of  the  propellant  particles,  the  buildup  of  a  strong  compaction  wave 
due  to  the  high  pressure  product  gases  and  eventual  shock  initiation  of  the  uoreacted 
explosive  upstream  of  the  ignition  front.  The  latter  part  of  th’s  proposed  scheme  of 
events  is  SDT  (shock-to-detonation- transition),  a  well-studied  method  of  initiating 
homogeneous  and  heterogeneous  explosives  without  the  use  of  a  thermal  ignition 
source.  In  SOT,  the  internal  energy  rise  accompanying  the  shock  front  is  Sufficient 
to  initiate  chemical  decomposition  of  the  unreacted  material. 

In  related  experimental  research,  Pilcher,  Beckstead,  Christensen  and  King  [A] 
investigated  the  effects  of  tube  deformation  on  OSDT  in  porous  beds  of  high  energy 
solid  rocket  propellants  (HMK-basedl.  The  research  involved  a  series  of  ten  OSOT 
experiments  in  addition  to  CO r r e s pond i n g  numerical  predictions.  HONDO  [5],  a  dynamic 
structural  analysis  code  provided  by  Sandia  Laboratories  was  used  *n  conjunction  with 
the  Krier/Van  Tassel  1  [6]  two-phase  flow  code  in  order  to  couple  flame  spreading 
through  the  porous  bed  with  the  motor  case  deformation.  For  the  material  examined 
( H M x ) .  confinement  of  the  product  gases  was  found  to  be  important  to  OSOT  outcome. 


Experiments  such  as  these  provide  useful  Information  about  the  stress  levels  and 
detonation  velocities  obtained  in  OSOT  events.  However,  in  order  to  fully  understand 
the  physical  mechanisms  behind  the  shock  formation  and  propagation  into  the  porous 
material,  one  cannot  rely  on  experiments  alone. 

The  work  to  be  presented  in  this  paper  represents  an  effort  to  analyze  and 
numerically  solve  a  mathematical  model  of  the  sequence  of  events  observed  in  the 
experiments  descirbed  above.  It  represents  an  important  extension  of  research  done 
earlier  by  Butler,  Krier  and  lembeck  [ 7  j  on  OSOT  in  HMX .  The  system  of  coupled 
partial  differential  equations  describing  the  one-dimensional ,  two-phase  reactive 
flow  will  be  presented  along  with  a  finite  difference  technique  used  to  numerically 
solve  them.  In  addition  to  this,  constitutive  equations  describing  the  states  of  the 
unreacted  solid  and  gaseous  products  are  also  given. 

Two  distinct  OSOT  mechanisms  will  be  discussed.  The  first  is  not  a  OSOT  in  the 
true  sense  of  the  word.  That  is,  a  steady  state  combustion  wave,  being  fueled  by  the 
rapid  burning  of  propellant  grains,  will  propagatae  through  the  granular  bed  and 
exhibit  CJ  detonation  properties.  However,  because  of  the  open  pores  and  rapid 
decomposition  of  the  solid  grain,  a  shock  wave  does  not  propagate  ahead  of  the 
detonation  front  as  found  in  detonation  waves  travelling  through  homogeneous 
materials.  The  second  OSOT  mechanism  does  involve  shock  heating,  but  only  after  the 
granular  bed  has  been  compacted  and  the  flow  of  hot  gases  has  been  terminated. 

THE  MODEL 

The  conservation  equations  (one-di mens i ona 1  !  describing  the  two-phase  (solid- 
gas)  flow  of  a  reactive  material  consist  of  a  system  of  six,  nonlinear,  coupled  par¬ 
tial  differential  equations. 


For  the  gas  phase  the  conservation  equations  for  mass,  momentum  and  energy  are 
expressed  in  the  Eulerian  form  as 
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and  for  the  solid  phase 
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The  terms  S,  M,  and  E  represent  the  interphase  mass,  momentum  and  energy  transfer 
terms  given  by 


gas  phase  transfer  terms 
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m  «  rup  -  C  (7b) 

£  *  r(Echeffl  -  u2p/2 )  -  Dup  -  Q  (7c) 

solid  phase  transfer  terms 

S  =  -r  (8a) 

R  *  C  -  r up  (8b) 

E  *  n.jp  *■  Q  -  r  u2p/2  (8c) 


Her®,  the  terms  EqT  and  EpT  represent  the  total  (sum  of  kinetic  and  internal)  ener¬ 
gies  in  the  gas  and  sol id^phases .  The  subscripts  g  and  p  denote  gas  and  particle, 
respectively.  In  Eqs.  (l)-(6),  the  phase  densities  and  pj  are  defined  as 

°1  *  °g*  and  °2  s  ^ 

The  variable  +  represents  the  percent  by  volume  occupied  by  the  gas  phase,  a  * 
Vg/Vtotal*  The  state  equation  for  the  gas  phase  is  a  nonideal  Abel  equation 
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where  i s  a  covolume  term.  Tor  the  solid  phase  we  used  a  form  of  the  Mi  e-Gr;inei  sen 
equat i on . 
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In  Equation  (10),  PH$  and  eHr,  the  Mijgoniot  pressure  and  Hugoniot  energy,  are  used  as 
the  reference  state  and  r(v?)  is  the  Gruneisen  coefficient  representing  the  thermo¬ 
dynamic  derivative, 


r  (vs)  =■  vs  ( a p / a e ) v 


(11) 


In  Eqs.  7a-7c  and  8a-8c  the  term  r  -epresents  the  rate  of  mass  generation  per 
unit  volume  per  unit  time  and  is  given  by 

r  1  ~  ( l  -  a  )o  r  .  (U) 

P  p 

Here,  r  is  the  instantaneous  oarticle  radius  and  r  is  the  surface  burning  rate  spe¬ 
cified  as  a  function  of  pressure. 
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NUMERICAL  SOLUTION  TECHNIQUE 


The  governing  equations  (Eqs.  1-6)  discussed  in  the  previous  section  are  hyper¬ 
bolic  and  highly  nonlinear  due  to  the  chemical  reaction  terms.  A  method  of  lines 
technique  [8]  was  implemented  to  solve  the  system  of  differential  equations  on  a  high 


speed  computer  (Cyber-175).  The  technique  requires  that  each  spatial  derivative  in 
the  equations  be  expressed  as  a  second  order  (or  higher)  finite  difference  analog  at 
each  time  level.  The  resulting  set  of  ordinary  differential  equations  (DDE's)  are 
then  marched  in  time  by  a  standardized  DOE  solver.  Because  the  resulting  equations 
are  very  stiff,  we  chose  to  solve  them  using  Gear's  bacxward  differentiation  formulae 
( BO F )  routine  [9], 

If  one  discretizes  the  propellant  bed  length  into  NNI  cells,  the  number  of  dif¬ 
ferential  equations  needed  to  be  solved  results  in  6  x  (NNI).  Each  of  Eqs.  1-6 
should  be  solved  at  each  of  the  NNI  grid  points  for  every  time  level.  However,  in  a 
one- d i mens i ona 1 ,  one-directional  wave  propagation  problem,  some  computational  short¬ 
cuts  can  be  taken.  Since  the  wave  motion  initiates  at  one  end,  it  is  not  necessary 
to  solve  the  equations  over  the  entire  spatial  domain.  In  order  to  minimize  the 
computational  time,  we  introduced  a  technique  which  computes  the  x-location  of  the 
lead  compression  wave  after  each  time  step.  Any  cell  upstream  of  this  point  has  yet 
to  be  disturbed  by  the  flow  field.  Therefore,  it  is  only  necessary  to  integrate  from 
x*0  up  to  and  including  the  location  of  the  lead  wave.  As  a  test,  we  solved  a  base¬ 
line  case  using  both  techniques.  The  first  case  solved  all  6  x  NNI  equations  and  tne 
second  solved  only  the  “dynamic"  domain.  A  RMS  error  evaluation  of  the  pressure 
profiles  at  t  »  10  usec  gives  an  prr or  of  RMS  *  3.21  between  the  two  cases.  The  test 
case  which  solved  the  entire  6  x  NNI  equations  took  19B  cpu  seconds  where  the  case 
which  solved  the  “dynamic"  domain  only  took  98  cpu  seconds,  a  reduction  of  50X  in 
computer  time. 

RESULTS 

As  a  baseline  test  case  we  chose  to  model  a  packed  bed  of  granulated  propellant, 
20  cm  in  length  and  closed  at  both  ends.  The  propellant/explosive  was  hmx  (octogen) 
and  was  assumed  to  be  consisting  of  unisized  spheres.  Nonspherical  particles  can 
easily  be  modeled  by  simply  adjusting  the  sphericity  factor,  x/r  in  the  decompo¬ 
sition  rate  (r  *  (x/rp)(p2  f)>.  for  spheres,  x  *  3.  v 

The  initial  conditions  for  this  particular  case  include:  porosity  *  *  0.30, 
propellant  particle  size  r  0  *  100  um,  material  density  (HMX)  o0  *  1.90  g/cc  and 
initial  temperature  T  »  30u  K.  The  linear  burning  rate  for  the  propellant  particles 
is  f  *  2.54  (.001  P  •’)  cm/s  where  P  is  the  surrounding  pressure  measured  in 

(psi).  In  order  to  start  the  transient  combustion  process,  the  material  located  at  x 
=  0  is  assumed  to  be  ignited  at  t  »  0.  This  initial  configuration  is  typical  of  DSDT 
experiments  run  in  the  laboratory  [2]. 

figures  2-5  illustrate  the  profiles  of  various  important  parameters  as  time 
progresses  and  the  combustion  wave  propagates  through  the  material.  Note  that  the 
pressure  (fig.  2),  velocity  (fig.  31,  and  mass  generation  rate  (fig.  4)  all  asymptote 
to  steady  state  values  after  approximately  80  usec.  This  occurs  when  the  ignition 
front  is  located  approximately  15-17  cm  from  the  ignition  source.  This  is  typical  of 
run-up  lengths  observed  in  actual  0S0T  experiments. 

Figure  2  illustrates  the  pressure  profiles  obtained  in  the  propellant  bed  at  2.5 
usec  increments.  Because  of  the  pressure  scale,  the  plots  for  0<t<20  usec  have  been 
excluded  from  this  figure.  The  value  of  peak  pressure  is  shown  to  approach  17  GPa  in 
the  steady  state  limit. 

Once  the  steady  state  detonation  wave  has  developed,  the  bed  can  be  described  as 
three  distinct  zones  as  is  evident  from  all  the  figures.  The  region  upstream  of  the 
shock  front  is  undisturbed  and  remains  at  the  initial  thermodynamic  state.  At  the 
shock / i gn i t i on  front  the  propellant  rapidly  decomposes  within  a  very  thin  zone  of 
reaction.  This  Is  evident  from  the  sharp  spikes  in  the  mass  generation  rate  plots 
(Fig.  4).  As  in  any  steady  state  detonation  wave,  the  rapid  reactant  *  product 
transformation  provides  the  necessary  driving  force  for  the  steady  state  wave, 
finally,  downstream  of  the  reaction  zone  Is  an  all-gas  expansion  wave,  extending  back 
to  a  stationary  zone  of  expanded  product  gases. 

, In  Figure  5  the  porosity,  a  (gas  volume/total  volume)  of  the  mixture  is  plotted 
as  a  function  of  x,  the  bed  location.  4s  the  reaction  front  travels  through  the  two- 
phase  mixture,  the  solid  reactants  decompose  and  thus  the  value  of  a  increases  until 
a  ♦  1  and  an  all-gas  composition  exists.  An  interesting  observation  made  from  the  a 
-  x  profiles  is  the  compaction  wave  developed  at  the  early  times,  45  <  t  <  65  usec. 

It  will  be  shown  later  how  this  decrease  in  a  can  effect  the  entire  DSDT  process. 

That  is,  under  certain  circumstances  the  value  of  a  may  approach  zero,  indicating  a 
complete  blockage  of  the  product  gas  flow. 
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Fig.  2  Fre$sure-x  profiles  showing 
shock  development  and  transi¬ 
tion  to  steady  state  detona¬ 
tion. 


Fig.  3  Gas  chase  velocity  profiles 
for  same  case  as  shown  in 
Fig.  2. 
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Fig.  4  nat.e  of  oas  generation  from 
solid  phase  to  product  gas 
pnase  (Fig.  2). 


Fig.  5  7olu~e  fraction  of  gas  pnase 
snowing  compaction  in  early 
ti-es  (50rt'65  _s( . 


figures  2-5  show  the  ignition  'root  to  accelerate  at  f irst  ard  then  to  stabilize 
at  a  detonation  velocity  Dr  j  *  5  . 7  !"n/jsec.  The  accompanying  r  o  apman-,j  ou  g  u  et  ores- 
Sure  is  Orj  ,  17.0  GPa.  These  values  are  close  in  agreement  with  one -d i me n s i on  a  1 , 
planan  detonation  theory  for  granular  explosives  and  with  values  obtained  from  TIGER 
predictions  for  HMX. 

TIGER  values  for  this  case  are 
?CJ  -  17.3  GPa 

and  0(-j  *  7,05  mm/  u  s  ec 

The  values  of  CJ  pressure  for  various  loading  densities  are  shown  in  f  \  g .  6t 
Here,  the  loading  density  is  defined  as  00  *  (1  -  *)o„  where  0-  is  t h e  crystalline 
density.  Data  for  this  plot  were  taxer  from  TIGER  predictions,  Johansson  [ 10  ] ,  and 
the  predictions  made  with  our  ende.  411  three  sources  show  detonation  pressure  to  be 
a  linear  function  of  loading  density  squared. 

Before  proceeding  with  the  discussion,  some  comments  should  he  made  regarding 
the  general  structure  of  the  steady  state  wave  solutions  generated  by  our  code. 

Unlike  the  classical  ZNO  detonation  wave,  our  results  do  not  show  shock  compression 


of  the  upstream  reactants.  In  our  results  the  reaction  front  and  compression  front 
appear  to  travel  as  one.  This  Is  due  to  the  two-phase  nature  of  the  problem,  when  a 
detonation  wave  travels  through  a  homogeneous  material  only  stress  waves  can  he  pro¬ 
pagated  upstream  of  the  reaction  front.  This,  of  course,  neglects  heat  transfer  by 
conduction  which  Is  of  a  lower  order  time  scale  than  the  detonation  velocity.  How¬ 
ever,  in  a  two-phase  mechanical  mixture  both  stress  waves  and  hot  product  gases  can 
propagate  through  the  upstream  material,  thus  changing  the  structure  of  the  wave,  and 
the  detonation  mechanisms. 

In  any  case,  the  expansion  wave  following  the  steady  state  detonation  should 
satisfy  certain  fundamental  relations  of  fluid  mechanics.  In  particular,  any  point 
along  the  expansion  wave  is  related  to  the  CJ  point  by  [11]. 


CJ 
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(v-l)(u-uCj)/2cCJ 


,rrrr 


(17) 


where  r  is  the  polytropic  exponent,  u  is  the  particle  velocity  and  c  is  the  sound 
velocity.  It  is  also  known  that  the  particle  velocity  must  equal  zero  at  the  fixed 
boundary  condition.  Thus,  the  steady  state  pressure  at  the  fixed  wall  should  satisfy 
Eq.  17  with  u  *  0.  For  the  case  shown  in  Figs,  2-5,  P  is  calculated  to  be  4.6 

GPa  (46  kbar).  This  is  consistent  with  the  value  predicted  by  the  code,  thus 
verifying  that  an  actual  steady  detonation  has  been  predicted. 

4  second  observation  made  of  the  predicted  results  and  consistent  with  theory 
is,  (dP /dt )  x _g-0  once  the  steady  wave  has  developed  downstream.  This  indicates  that 
the  expansion  wave  reaches  a  zero  velocity  point  somewhere  between  the  detonation 
front  and  rear  boundary  and  thus,  the  pressure  at  the  rear  boundary  remains  constant 
with  time.  It  has  been  shown  in  the  HOT  calculations  given  in  Ref.  7  that  in  some 
cases  (moderate  burning  rate  propellants)  high  pressures  can  be  predicted,  but  a 
steady  state  detonation  does  not  develop.  (This  is  still  a  hazard  in  a  solid  rocket 
motor  since  the  pressures  attained  well  exceed  the  strength  of  the  motor  casing.) 
However,  if  the  porous  b»d  was  ignited  with  the  purpose  of  developing  a  certain 
strength  shock,  say  to  shock  initiate  a  downstream  solid  charge  [14],  then  the  DSDT 
has  failed,  unless  Equation  (17)  is  satisfied. 

PLUG  FORMATION 

In  some  instances  the  porosity  upstream  of  the  convective  front  may  approach 
zero  (1001  solid).  Our  results  show  this  to  occur  early  in  the  process  before  the 
detonation  wave  develops.  Also,  it  is  more  likely  to  occur  when  the  propellant 
particles  are  large  and  slow  burning  initially.  When  this  occurs  the  flow  of  hot 
product  gases  is  restricted  and  thus,  the  upstream  material  can  no  longer  be  heated 
by  convective  means.  The  pore  collapse  and  plug  formation  are  illustrated  in  Fig. 

7.  Part  A  represents  the  initial  con f  i  gu ra t i on  ,  just  after  ignition,  and  Part  R 
illustrates  the  plug  formation.  Once  the  plug  has  formed,  the  material  downstream  of 
the  plug  is  still  reacting,  acting  as  a  driving  force  for  the  upstream  matrix  of 
material.  This  c on f i gu r a t i on  is  similar  to  DOT  experiments  performed  on  cast 
explosives  [12]  where  one  end  is  confined  and  end-burning  initially.  The  result  is 
coalescence  of  stress  waves  upstream  of  the  reacting  zone  and  eventual  shock 
initiation. 


Figure  I  is  a  plot  of  predicted  a  -  x  profiles  for  a  packed  bed  of  hmx.  The 
particles  are  spheres  with  a  200  um  diameter  and  the  initial  porosity  is,  a  *  0.32. 
The  profiles  are  given  at  2.5  usee  time  increments.  As  shown  in  the  figure,  the 
porosity  approaches  a  *  0.10  at  t  *  75  usec.  At  that  time  the  plug  was  assumed 
formed. 

The  sequence  of  events  occurring  after,  a  plug  has  formed  are  shown  for  a 
different  case  in  Fig.  9.  Here,  the  str’ss  -  x  profiles  are  plotted  as  a  function  of 
t  ,  the  time  after  plug  formation.  The  compression  wave  generated  from  the  reacting 
grains  is  reflected  off  of  the  plug  and  a  strong  shock  wave  is  transmitted  to  the 
upstream  material.  Since  the  code  does  not  yet  have  the  capability  to  model  the 
thermal  explosion  of  a  reactive  material  which  has  been  shocked,  one  can  only 
estimate  a  detonation  run-up  length  from  experimental  Pop*plot  data  on  similar 
materials.  For  example,  for  a  shock  wave  of  more  than  1.0  OP  a  Dick  [13]  has  shown 


the  run-up  length  to  be  less  than  0.5  cm  for  1.24  g/cc  HMX  (* 
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CONCLUSIONS 

The  work  discussed  here  addresses  the  problem  of  reactive,  two-phase  gas 
dynamics  associated  with  the  development  of  a  detonation.  It  has  been  pointed  out 
that  as  new,  more  energetic  solid  propellants  are  proposed  for  new  propulsion  systems 
the  ability  to  understand  the  physics  of  this  hazard  become  more  necessary. 

The  research  of  flow  simulation  through  modeling  can  give  (and  has  provided) 
Information  that  more  clearly  explains  the  transient  events.  Or,  put  more  strongly, 
It  would  be  virtually  meaningless  to  attempt  to  understand  DSOT  by  performing 
experiments  alone,  because  it  would  he  impossible  to  ascertain  the  flow  coupling 
which  provides  the  observed  pressures  and  flame  front. 

We  have  shown  two  distinct  ways  in  which  a  bed  of  granulated  high  energy 
propellant  can  detonate  from  a  controlled  deflagration.  In  both  cases  the  bed  is 
initially  confined  and  has  a  high  solids  loading. 

At  first  we  showed  that  a  bed  of  granular  particles,  very  small  n  size  (100 
um),  can  exhibit  detonation  characteristics  by  burning  in  a  convective  manner.  Here, 
the  important  process  was  rapid  gas  generation  from  the  burning  particle.  The 
detonation  characteristics  were  also  shown  to  match  very  closely  those  predicted  by 
the  TIGER  code  and  experiment. 

The  second  part  of  our  proposed  research  in  OSDT  showed  that  if  a  plug  (»  ♦  0) 
forms  in  the  granular  bed,  prior  to  the  rapid  convective  acceleration,  the  upstream 
plugged  material  can  be  shock  loaded  to  strengths  great  enough  to  initiate  detonation 
within  a  few  millimeters. 

In  summary,  the  research  shows  a  unique  complex  inter-coupling  between 
propellant  properties,  (such  as  density,  porosity,  energy  content,  burning  rate)  and 
flow  containment  (such  as  permeability,  convective  heat  transfer,  and  grain 
deformation)  that  occur  in  a  0S0T  event.  A  recent  report  by  Baer  and  Nunziato  [1*] 
confirms  that  DSOT  modeling  of  the  type  described  herein  can  accurately  model 
detonation  experiments  in  porous  reac t i ve  so  1 i ds . 
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fig.  7  Illustration  of  second  branch  of  DSOT. 

Part  A  represents  initial  deflagration 
state.  Part  B  shows  the  plug  formation 
and  upstream  compaction. 
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Fig.  6  Detonation  pressure 
as  a  function  of 
initial  loading  den¬ 
sity  squared  for  HMX . 


.  loooe-n 


Fig.  8  Volume  fraction  of  gas  phase  showino 
complete  plug  formation  ( t>— 0)  for  HMX  , 
200  urn. 


Fig.  9  Shock  propagation  upstream  of 
plug  after  plug  formation. 

The  plug  is  formed  around  x  *  10  cm. 
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TABLE  OF  INPUT  PARAMETERS 

b  0.001  in/s  (psi)n 

^ chem  6-2 

fpg  fpq  •  ( ( l-*)/e)2(276  ♦  5(Re/(l-*))-87  (Ref.  7) 


cvg 

1.57  x  107 

erg/g-K 

cvp 

1.00  x  107 

erg/g-K 

n 

0.90 

Pr 

0.7? 

ro 

100  -  200 

I*  Ifl 

r0 

1  .65 

1 

n 

4.00  cc/g 

*0 

0.30 

u  9 

1.8  x  10'4 

g/cm-s 

L 

20  cm 

Po 

1.0  x  106 

dynes/cm 

°Po 

1.9  g/cc 

N  N 1 

100  cells 

To 

300  K 

aEign 

5.0  x  107 

erg/g 
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TRANSITION  TO  DETONATION  FROM  RAPI0-CDMP9ESS10H 
(RAMP-WAVES)  GENERATED  IN  A  BURNING  POROUS  BED 

*  «* 

Christopher  A.  Cudak  ,  Herman  Krier  and  P.  Barry  Rutler' 

University  of  Illinois  at  UrOana-Champaign 
tlrbana,  Illinois 

ABSTRACT 

Increasing  the  nitramine  content  of  solid  rocket  propellants  Increases  the  overall  performance 
of  the  system  as  well  as  the  sensitivity  to  detonation  by  shock  initiation.  In  some  Instances  a 
confined  zone  of  granulated  propellant  adjacent  to  a  zone  of  cast  propellant  can  provide  a  rapid 
enough  pressure-rise  rate  to  shock  initiate  the  cast  material.  If  the  cast  propellant  is  porous, 
the  detonation  will  initiate  at  some  location  ahead  of  the  granulated  bed/cast  material  Interface. 
The  work  presented  here  is  an  effort  to  numerically  model  this  Deflagration  to  Shock  to  Detonation 
Transition  (OSDT)  event.  Results  will  be  presented  showing  the  detonation  build  up  for  propellant 
beds  with  various  initial  configurations  and  boundary  conditions. 

LIST  OF  SYMBOLS 

Cv  specific  heat  (erg/g/K) 

e  specific  internal  energy  (erg/g) 

E,  total  energy  (erg) 

E  activation  energy  (erg/mole) 

G  shear  modulus  (dynes/cm2) 

P#  pressure  (dynes/cmz) 

P  maximum  axial  stress  (dy/cm*) 

0  chemical  energy  release  rate  (erg/g/s) 

R  gas  constant  (erg/g/K) 

t#  time  (s) 

t  characteristic  rise  time  (s) 

T#  temperature  (K) 

T  characteristic  burn  temperature  (K) 

u  particle  velocity  (cm/s! 

v  specific  volume  (cnr/g) 

V  volume  (cm3) 

U  decomposed  mass  fraction 

x  spatial  location  (cm) 

Y  yield  strength  (dynes/cmz) 

l  frequency  factor  (1/s) 


INTRODUCTION 

The  probability  of  a  detonation  occurring  in  a  solid  propellant  rocket  motor  greatly  increases 
when  secondary  high-explosives  are  used  as  constituents  in  the  propellant  mixture.  High-energy 
nitramines  such  as  HMX  (octogen)  are  often  included  in  the  propellant  formulation  in  order  to 
Increase  the  specific  impulse  and  thus  overall  performance  of  the  rocket  motor.  However,  when 
utilizing  these  high-energy  formulations  the  hazard  of  a  Deflagration  to  Shock  to  Detonation 
Transition  (DSOT)  becomes  more  likely.  A  OSDT  event  is  defined  as  a  controlled  (subsonic) 
accelerating  deflagration  wave  making  a  transition  to  a  high  order  steady  detonation  wave.  Once  a 
region  of  granulated  material  has  been  Ignited  in  a  confined  configuration,  the  OSDT  can  occur  in 
the  order  of  tens  of  microseconds  and  the  result  is  total  destruction  of  the  rocket  motor 
assembly. 

Although  little  data  is  available  on  the  actual  OSOT  process,  most  researchers  agree  that  in 
order  for  OSDT  to  occur  the  rocket  motor  grain  must  first  be  damaged.  Gram  fragmentation  can  be 
the  result  of  a  handling  accident,  or  case  and/or  nozzle  failure  during  operation.  The  high 
surface-to- volume  ratio  fragments  provide  for  the  rapid  pressurization  rate  necessary  to  shock 
initiate  the  remaining  grain.  Thus,  by  increasing  the  nitramine  content  of  the  propellant  mixture 
the  propensity  to  detonate  is  increased  in  two  ways;  the  decomposition  rate  of  the  propellant 
Increases  and  the  shock  sensitivity  Increases. 
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In  this  pap er  we  propose  two  possible  methods  In  which  OSOT  can  occur  In  a  solid  rocket 
motor.  In  the  first  case  a  section  of  the  motor  grain  Is  completely  fragmented,  and  In  the  second 
case,  the  grain  Is  only  partially  fragmented.  OSOT  In  a  fully  granulated  bed  (Case  1)  was  the  topic 
of  several  research  papers  Cl,2]  by  our  group  and  will  be  discussed  here  only  to  provide  a 
comparison  with  the  second  scenerlo  proposed.  In  Casa  1,  the  transition  to  detonation  occurs  within 
the  granulated  material  whereas  In  Case  2  It  occurs  In  the  upstream  cast  material.  For  the  second 
scenerlo  the  granulated  zone  Is  not  long  enough  to  OSOT,  however,  the  nonfragmented  material 
upstream  can  still  detonate  by  shock  initiation.  In  both  cases,  the  confined  burning  of  the 
propel lant/exploslve  fragments  Is  irftat  drives  the  deflagration  wave  to  detonation.  The  reader 
should  consult  Ref.  3  for  a  general  overview  of  OSOT  in  different  propellants  and  explosives  and 
Refs.  4*5  for  experimental  results.  Before  presenting  the  analysis  and  modeling  of  the  second  case, 
a  brief  review  of  Case  1  Is  appropriate. 

CASE  1 

Teg In  by  considering  the  rocket  motor  shown  In  Fig,  1.  For  Illustration  purposes  we  have 
selected  a  center-burning  grain  configuration.  As  far  as  we  know,  the  propellant  will  not  detonate 
unless  a  crack,  filled  with  broken  fragments.  Is  developed  in  the  cast  grain  [3].  That  Is  to  imply, 
the  rapid  pressurization  necessary  to  shock  Initiate  the  cast  propellant  cannot  be  generated  by 
extra  "side-wall*  surface  burning  (conductive  combustion)  alone.  Therefore,  to  Illustrate  a  OSOT 
event  we  have  assumed  there  Is  a  crack  In  the  propellant  and  that  within  the  crack  we  have  assumed  a 
packed  bed  of  propellant  fragments  of  known  surface-to- volume  ratio  to  exist.  The  fragment  filled 
crack  Is  enlarged  and  shown  In  Fig.  2.  At  some  arbitrary  Initial  time  t0,  the  bed  Is  assumed 
burning  at  x  •  0  (left-end  in  figure).  Because  of  the  assumed  high  surface-to- volume  ratio  of  the 
confined  sub-millimeter  size  particles,  product  gas  generation  Increases  beyond  the  level  necessary 
for  steady-state  motor  operation  indicating  that  the  amount  of  gas  being  generated  by  the 
decomposing  propellant  far  exceeds  the  amount  exiting  the  nozzle.  Under  these  conditions,  the  crack 
in  the  motor  can  be  modeled  as  a  bed  of  granulated  propellant  confined  on  both  ends  (Fig.  2).  As 
time  progresses,  pressure  gradients  are  developed  as  a  result  of  localized  burning.  Conserving 
momentum,  the  pressure  gradients  act  to  drive  the  hot  product  gases  Into  the  cracks  and  exposed 
ptjres  between  the  propellant  fragments.  Convective  heat  transfer  from  the  hot  gas  to  the  propellant 
surface  will  then  Ignite  more  of  the  particles,  eventually  leading  to  shock  compression  of  the 
upstream  propellant  and  subsequent  detonation.  Figure  3  shows  a  steady-state  pressure  profile  some 
time  after  detonation.  Ahead  of  the  detonation  wave  are  unreacted  fragments  of  propellant.  Within 
the  detonation  front  the  grains  are  rapidly  consumed  and  to  the  downstream  side  of  the  wave 
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Fig.  3.  Pressure-x  profile  of  steady  state 
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Fig .  4.  Sketch  of  granulated  bed/cast  bed 
configuration  (Case  2)  for  OSOT. 
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Fig.  6.  Pressure-rise  rates  in  granular  bed 
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the  product  gases  expand  in  a  Taylor  expansion  wave.  This  process  is  referred  to  by  some  as  the 
accelerated  convective  bum  model  and  is  described  in  detail  in  the  model  developed  by  Butler, 
Lembeck  and  Krier  [1]. 


CASE  2 

5”  second  OSOT  scenerlo  Involves  a  region  of  granulated  propellant  providing  the  Impetus  to 
shock  initiate  an  adjacent  region  of  cast  explosive.  This  is  illustrated  in  Fig.  4.  The  cast 
material  (Zone  1)  can  contain  'blind*  pores,  but  is  assumed  to  be  Impermeable  to  the  flow  of  hot 
gases  from  the  granular  zone  (Zone  2).  This  implies  that,  unlike  the  first  OSOT  scenerlo  discussed, 
only  stress  waves  can  be  transmitted  across  the  Zone  1/Zone  2  Interface.  Figure  5  shows  a  schematic 
representation  of  the  sequence  of  events  leading  to  this  type  of  DSDT.  Superimposed  on  each  section 
of  the  figure  is  a  solid  line  representing  the  local  gas  porosity  (volume  of  gas/total  volume)  as  a 
function  of  x,  the  bed  location.  A  value  of  »  equal  to  unity  represents  a  zone  of  all  gas  while  j 
equal  to  zero  indicates  a  homogeneous  solid. 

Part  A  shows  a  burning  granular  zone  adjacent  to  a  porous,  east  explosive.  Here,  the  heavy 
black  dots  are  representative  of  microvoids  in  the  cast  material.  Illustrated  in  Part  B  is  the 
collapse  of  the  pores,  a  result  of  the  stress  load  transmitted  across  the  granulated  bed/cast 
explosive  Interface.  Parts  C  and  0  show  the  length  of  the  pore  collapse  zone  to  increase  with  time 
as  the  lead  compression  waves  travel  farther  into  the  explosive.  The  finite  compression  waves 
coalesce  into  a  shock  front  which  then  shock  initiates  the  cast  explosive  downstream  of  the 
Interface.  From  this  location  a  detonation  wave  propagates  through  the  porous  material  while  a 
retonation  wave  propagates  back  through  the  compressed  material  (Part  E). 

In  related  research,  Setchell  [6]  studied  ramp-wave  initiation  mechanisms  using  waves  with 
rise-times  of  0.3  and  0.8  us.  In  Ref.  6  the  ramp-waves  were  developed  by  propagation  of  an  impact- 
generated  shock  wave  through  a  material  which  had  the  unique  property  of  spreading  out  the  shock 
front  into  a  ramp-wave.  In  the  work  presented  here,  the  rise-times  were  much  longer,  typically  of 
the  order  of  tens  of  microseconds.  These  slower  rise-times  are  typical  pressure-time  histories 
generated  within  an  accelerating  deflagrating  porous  bed. 
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ANALYSIS  OF  OSOT  IN  A  PARTIALLY  GRANULATED  AED 

In  a  recent  AGARD  paper  [2]  we  showed  that  under  certain  conditions,  DSDT  can  occur  in  a 
granulated,  porous  bed  of  HMX  (Case  1  scenerio),  A  few  of  the  necessary  conditions  for  DSOT  to 
occur  identified  in  Ref.  2  are:  sufficient  product  gas  confinement,  a  high  initial  loading  density, 
a  rapid  gas  generation  rate,  a  bed  length  greater  than  the  critical  detonation  run-up  length  and 
submit limeter  sized  particles.  The  steady  state  detonation  solutions  predicted  in  Ref.  2  were 
obtained  by  numerically  solving  the  defining  system  of  time-dependent  conservation  equations  in 
conjunction  with  the  proper  constitutive  relations.  Values  predicted  for  CJ  pressure,  CJ 
temperature,  detonation  velocity  and  detonation  run-up  length  were  in  close  agreement  with  values 
predicted  by  the  thermoequilibrium  code  TIGER  [7]  and  experimental  data.  The  reader  is  referred  to 
Ref.  2  for  more  on  this  subject. 

The  scope  of  the  work  presented  in  this  paper  is  to  model  DSDT  in  the  granulated  bed/cast 
explosive  configuration  shown  in  Fig.  4.  It  is  assumed  that  the  granulated  bed  is  not  long  enough 
to  undergo  DSDT,  but  by  convective  burning  provides  the  driving  force  necessary  to  shock  initiate 
the  upstream  cast  explosive.  For  this  analysis  the  cast  explosive  was  assumed  to  have  'blind'  pores 
and  to  be  impermeable  to  the  flow  of  product  gases  generated  in  the  granular  bed.  This  implies  that 
only  stress  waves  can  be  transmitted  across  the  granular  bed/east  explosive  Interface. 

The  first  step  to  modeling  OSOT  in  the  Case  2  configuration  is  to  determine  the  rate  at  which 
the  reacting  granulated  bed  stresses  the  cast  explosive.  This  task  was  accomplished  by  running  the 
OSOT  code  discussed  in  Ref.  2  for  granular  beds  with  lengths  less  than  the  detonation  run-up  length 
and  recording  the  pressure-rise  rate  at  the  end  opposite  the  igniter.  Plotted  in  Fig.  6  are  the 
P  -  t  functions  for  propellant  beds  of  various  particle  sizes.  As  shown  in  the  figure,  the 
pressurization  rate  in  the  granular  bed,  dP/dt,  is  strongly  dependent  on  the  size  of  the  particles 
being  consumed.  The  larger  the  particles  being  consumed,  the  slower  the  pressurization  rate  on  the 
interface.  Thus,  from  this  type  of  analysis  we  have  been  able  to  determine  the  P(t)  function  used 
as  a  boundary  condition  for  the  stress  wave  analysis  in  the  cast  explosive  upstream. 

The  P  -  t  functions  obtained  from  this  modeling  effort  have  been  linearized  and  are  expressed 
as 

P(t)  -  (P*  -  P0)(t/t#)  ♦  P0  t  £  t*  (la) 

P(t)  •  P*  t  >_t*  (lb) 

where  P*  is  the  maximum  stress  in  the  bed.  The  parameter  t*  is  a  characteristic  rise-time  for  the 
ramp-wave  input  function. 


GOVERNING  EQUATIONS 

The  Lagrangian  or  material  form  of  the  governing  equations  are  incorporated  in  the  hydrodynamic 
analysis  of  a  continuous  material  with  a  moving  boundary.  In  the  problem  addressed,  the  moving 
ooundary  is  a  result  of  the  applied  stress  load  from  the  burning  granulated  bed.  The  one- 
dimensional  unsteady  form  of  the  conservation  of  mass,  momentum,  and  energy  equations  are  expressed 
for  the  total  mechanical  mixture,  as 


v  ,Jx 

(2) 

PX 

(3) 

-  P  vt  ♦  Q 

(4) 

In  the  above  expressions  v  represents  specific  volume,  u,  particle  velocity,  e,  the  specific 
internal  energy,  P,  the  total  stress,  and  0  the  heat  added  per  unit  mass  per  unit  time.  The 
subscripts  x  and  t  indicate  partial  derivatives  with  respect  to  the  Lagrangian  spatial  coordinate 
and  time  respectively.  In  addition  to  the  conservation  equations,  a  material  equation  of  state, 

P  ■  P(v,e),  is  needed  in  order  to  provide  for  mathematical  closure. 

For  the  solid  material  the  equation  of  state  P.(vs,e.)  and  the  caloric  equation  es(vs,Ts)  are 
expressed  in  terms  of  a  Helmholtz  free  energy  function  [8j,  u(vs,Ts)  and  its  thermodynamic 
derivatives  through  the  Second  Law  of  Thermodynamics  reciprocity  relations 


With  the  assumption  that  the  Gruneisen  coefficient  is  constant,  the  Helmholtz  free  energy  function 
takes  the  following  form  [8] 
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where  ?  Is  the  Gruneisen  coefficient  defined  by  the  thermodynamic  derivative 

r  ( v )  *  -  v  .£ 

'  ;ev 

and  Cvs  is  the  specific  heat  at  constant  volume  of  the  solid  phase 
C 


VS 


'^e^ 

TTv 


'71 


[i) 


(9) 


The  term  J(v  )  in  Eq.  7  is  a  nonlinear  volume-dependent  function  determined  from  shock  Hugoniot 
experiments  [8], 


With  the  introduction  of  product  gases  into  the  system,  a  constitutive  law  for  the  gas  phase 
must  also  be  provided.  A  nonideal  covolume  equation  of  state  was  chosen 

Pg  *  W1  *  3V  (10) 

where  R  is  the  gas  constant  and  3  is  a  covolume  correction  term.  The  value  of  3  is  determined  from 
the  values  for  pressure,  temperature  and  density  at  the  CJ  state  predicted  by  the  TIGER  code.  Table 
I  gives  a  listing  of  those  values  for  several  loading  densities  of  HMX.  In  accordance  with  the 
reciprocity  relations  defined  earlier  (Eqs.  5-6)  the  caloric  equation  of  state  for  the  gas  phase  is 

e9  *  *"vg  '^9  *  T90* 

where  Cvg  is  the  specific  heat  at  constant  volume  of  the  product  gases. 

REACTIVE  POROUS  MATERIALS 

The  governing  equations  for  conservation  of  mass,  momentum  and  energy  defined  in  the  previous 
section  (Eqs.  2-4)  are  expressed  in  terms  of  the  thermodynamic  properties  (P,v,e)  of  the  two-phase 
mixture  as  well  as  the  dynamic  variable  u.  Additional  relations  are  needed  to  separate  the 
individual  phase  properties  from  those  of  the  mechanical  mixture.  It  is  assumed  that  any  arbitrary 
volume  VT  within  the  continuum  can  contain  both  solid  and  gas  phases.  Thus,  the  individual  phase 
volumes  sum  to  equal  the  total 

VT  »  Vg  «•  Vs  or,  v;=  (l-w)v  +  Wv$  (12a) 

and,  likewise  for  total  energy  ET 

Et  *  Eg  *  Es  or,  eT= ( l-W)eg  +  We$  (12b) 

It  is  also  assumed  that  P  *  P$  *  Pg. 

We  define  the  material  porosity  as  the  ratio  of  total  volume  to  volume  occupied  by  the  solid 

phase 

i  *  VT/VS  (13) 


Initially,  the  volume  not  occupied  by  the  solid  material  is  assumed  to  be  "assies*.  Ay  introducing 
porosity  into  the  material,  the  mathematical  description  is  not  complete  unless  an  additional 
Independent  equation  is  also  provided. 


Carroll  and  Holt  [9]  have  performed  extensive  research  in  the  area  of  mathematically  modeling 
the  collapse  of  a  porous  material  under  an  applied  external  load.  In  their  model  the  porous  matrix 
is  treated  as  a  hollow  sphere  where  the  inner  and  outer  radii  are  chosen  such  that  pore  size  and 
overall  porosity  of  the  material  are  accurately  represented.  The  assumed  pore  collapse  occurs  in 
three  phases:  (l)  elastic  phase,  where  the  solid  elastically  deforms:  (2)  elastic-plastic  phase, 
where  plastic  deformation  begins  at  the  Inner  radius  and  progresses  outward  until  plastic 
deformation  begins  at  the  outer  radius:  and  (3)  plastic  phase  where  plastic  deformation  occurs 
throughout.  The  P  -  i  relations  for  the  three  phases  of  compaction  and  the  appropriate  range  over 
which  each  applies  are  given  by. 


elastic  phase  i  1  *  *1 
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where  the  limits  between  the  three  phases  are  given  by 
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with  Y  and  5  the  yield  stress  and  shear  modulus  respectively.  During  pore-collapse,  P  *  Ps. 

T*E  C0M8UST10H  MQ0EL 

A  first  order  Arrhenius  burn  model  is  used  to  describe  the  chemical  decomposition  of  the 
reactive  material.  It  is  given  as 

—  «  -  ZU  exp  -  £  /ITT  '  wnere,  14  ■  m^  /  mjo  (19) 

In  the  above  equation  u  is  the  mass  fraction  of  unreacted  explosive,  Z  is  the  frequency  factor,  E* 
is  the  activation  energy,  R  the  universal  gas  constant  and  T  the  characteristic  burn 
temperature.  During  compression  of  the  porous  bed,  T  represents  a  “hot  spot"  temperature  due  to 
irreversible  heating  at  pore  sites.  The  localized  hot  spot  temperature  is  different  from  the  bulk 
shock  temperatures  and  is  calculated  from  a  theory  developed  by  Hayes  [10]. 


The  underlying  assumption  in  this  theory  is  that  the  shocked  porous  material  is  at  one  of  two 
possible  temperatures,  a  bulk  shock  temperature  T,  or  a  hot  spot  temperature  TH.  The  energy 
deposited  by  the  shock  wave  is  equated  on  a  mass  fraction  basis  to  the  sum  of  the  reversible  work 
done  in  isentropical ly  compressing  the  bulk  of  the  material,  plus  the  irreversible  heating  of 
localized  hot  spots 
P  ♦  ? 

— (vT(j  -  v)  ■  WH  e(v,TH)  -  (1  -  WH)  e,  (P).  (20) 

In  Eq.  20  the  left-hand  side  represents  the  total  energy  deposited  in  the  material  by  the  shock  of 
strength  P.  The  term  e^P)  represents  the  energy  required  to  isentropical  ly  compress  the  bulk  of 
the  material  to  the  final  shock  pressure  and  the  remaining  energy  term,  e(v,TH),  is  the  energy 
available  to  irreversibly  heat  the  hot  spots.  The  Hayes  model  assumes  the  mass  fraction  of  the  hot 
spots,  Wh,  to  be  equal  to  the  preshock  volume  fraction  of  pores. 


WH  *  vTo/VSO  ■  1  *  -  1  (21) 

uere,  the  subscript  ‘To’  represents  the  initial  porous  state  and  the  subscript  'so1  refers  to  the 
rd-ogeneous  initial  state. 
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Pig.  7.  A  progress  diagram  for  a  stress 
wave  propagating  through  a  porous 
bed  of  rMX  (c-  *1.3,  .*1.4615), 
translating  t6°a  steady  detonation 
and  retonation  (P**3  5?a  and 
t**10  .sec) . 
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Fig.  3.  3hysical  plane,  showing  locus  com- 
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NUMERICAL  SOLUTION  TECHNIQUE 

Equations  2-19  completely  define  the  fluid  motion  anq  thermodynamic  state  of  a  continuous, 
porous  media.  The  system  of  equations  were  solved  by  using  a  finite  difference  numerical  technique 

patterned  after  the  WONOT  code  [ll].  At  t  •  0  the  bed  of  porous  propel!ant/e»plosive  is  discretized  .  9 

into  J  cells  labeled  from  left  to  right  as  j  *  l,  2,.  .  .  J.  The  thermodynamic  properties  pressure, 

temperature,  internal  energy  and  specific  volume  are  assumed  to  be  constant  over  the  width  of  eacn 

cell.  At  the  boundaries  between  the  cells  values  for  particle  velocity  and  spatial  location  are 

assumed  known.  The  reader  Is  referred  to  Ref.  12  for  a  listing  of  the  finite  difference 

approximations  to  the  governing  differential  equations  and  constitutive  relations. 

RESULTS 

As  stated  earlier,  the  purpose  of  this  work  Is  to  predict  the  transient  events  leading  to 
detonation  for  a  bed  of  porous,  cast  hmx  being  stress  loaded  at  one  end  by  a  burning  granulated  bed 
of  the  same  explosive.  In  our  preliminary  calculations  we  found  that  rise-times  of  10  <  t  <  30 

us  were  typical  of  most  of  the  granulated  beds  studied. 
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Figure  7  shows  the  stress-x  profiles  predicted  at  various  times  for  a  porous  bed  of  hmx.  This 
bed  is  ramp  loaded  at  Jhe  left  boundary  to  a  maximum  stress  of  3.0  GPa  and  the  characteristic  rise¬ 
time  in  tnis  case  is  t  «  10  u$.  As  the  left  boundary  is  toaded,  stress  waves  propagate  through  the 
porous  material,  collapsing  the  voids  and  coalescing  into  a  shock  front  at  approximately  x  «  5  cm. 

At  the  time,  t  *  24.7  ys,  the  irreversible  heating  of  localized  areas  by  the  shock  wave  causes  an 
explosion  to  occur.  The  rapid  release  of  reactant  chemical  energy  supports  the  detonation  wave 
which  propagates  in  the  +x  direction  at  a  velocity  of  7.1  mm/us.  At  the  same  time,  a  backward 
traveling  retonation  propagates  through  the  compressed  material  that  remains  between  the  detonation 
front  and  x  *  0  boundary.  This  wave  travels  at  a  slightly  greater  velocity,  3.3  mm/us  than  the 
detonation.  This  is  expected  since  the  material  through  which  it  is  propagating  is  precompressed 
and  therefore  has  a  higher  detonation  velocity.  The  values  predicted  by  our  code  for  detonation 
velocity,  0  *  7.1  mm/us,  CJ  pressure,  Prj  «  20.0  GPa  and  CJ  temperature,  TCJ  ■  450D#K  are  in  close 
agreement  with  values  predicted  by  the  tnermochemical  code  TIGER  (see  Table  I). 

Moreover,  the  transition  to  detonation  is  delineated  in  Fig.  8.  Here,  the  dashed  line 
represents  the  stress  wave  propagating  through  the  HMX,  initiating  the  detonation  at  x  *  4.37  cm. 

The  solid  lines  represent  the  detonation  front  locus  and  retonation  front  locus.  The  discontinuity 
in  slope  at  t  >  24.7  us  indicates  the  time  at  which  the  tranistion  occurs  and  the  corresponding 
length  is  termed  the  detonation  run-up  length,  i-j.  The  discontinuity  in  the  retonation  front 
profile  occurs  when  the  wave  Interacts  with  the  driving  piston. 

Experimental  data  are  not  available  for  10-30  us  rise-time  ramp-waves.  Therefore,  in  order  to 
validate  our  calculations,  we  have  extrapolated  our  data  in  the  10-30  vs  range  back  to  a  regime 
where  data  are  available.  Recall,  Setchell  [6]  did  experiments  with  t  »  0.3  us  and  Dick  [13]  has 
made  run-up  measurements  for  shock-to-detonation  transition  (SOT)  in  porous  HljJX  (OTo  •  1.24  g/cc). 

By  extrapolating  our  t  •  trj  data  for  oTo  •  1.30  g/cc  HMX  (Fig.  9)  down  to  t  •  0  us,  a  run-up 
length  of  -  5mm  is  predicted.  Extrapolating  the  Pop-plot  data  in  Ref.  13  gives  a  run-up  length 
of  —  3  . 

Finally,  a  comparitive  study  was  made  for  various  ramp-wave  rise-times  and  the  results  are 
shown  in  Fig.  10.  Here,  the  material  being  stress  loaded  has  the  same  initial  conditions  as  the 
case  studied  in  Fig.  7.  The  Important  point  here  Is  to  compare  the  detonation  run-up  length,  trg, 
for  the  various  ramp  waves.  The  longer  the  rise-time  is,  the  greater  is  the  run-up  distance.  This 
is  an  important  point  in  rocket  motor  hazard  analysis.  The  granulated  bed/cast  explosive 
configuration  is  representative  of  a  rocket  motor  which  has  been  partially  fragmented.  Even  though 
the  length  of  the  fragmented  propellant  is  not  long  enough  to  dSdT,  the  rapid  rise-rate  may  shock 
initiate  the  Intact  propellant.  However,  if  the  distance  needed  to  shock  initiate  the  cast  material 
is  greater  than  any  dimension  of  the  rocket  motor,  DS0T  is  Impossible.  The  motor  will 
overpressurize  and  rupture,  but  it  cannot  detonate. 

conclusions 

It  has  been  shown  by  numerical  simulation  that  under  certain  circumstances  a  porous,  cast 
explosive  will  transition  to  a  detonation  when  stress  loaded  in  a  ramp-wave  fashion  by  an  adjacent 
bed  of  granulated  explosive.  Rise-times  for  the  waves  studied  ranged  from  10  uS  to  30  us.  In  most 
cases,  the  bed  detonated  at  some  location  downstream  of  the  stress-loaded  boundary  and  an  energy 
sustaining  retonation  wave  was  shown  traveling  back  towards  the  driving  wall.  The  values  predicted 
for  the  detonation  properties  were  in  agreement  with  TIGER  data  although  the  run-up  lengths  were 
much  greater  than  any  predicted  for  SOT  experiments. 

At  the  present  time,  the  analysis  is  being  applied  to  propellant  beds  with  initial  porosities, 
1.05  <  x„  <  1.5.  In  addition,  the  sensitivity  of  the  ignition  criterion  to  0S0T  is  being  studied 
along  with  alternate  chemical  decomposition  formulations  (Eq.  19). 
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